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ABSTRACT 

Cold gas streaming along the dark-matter filaments of the cosmic web is predicted to 
be the major provider of resources for disc buildup, violent disk instability and star 
formation in massive galaxies in the early universe. We study to what extent these cold 
streams are traceable in the extended circum-galactic environment of galaxies via Lya 
absorption and selected low ionisation metal absorption lines. We model the expected 
absorption signatures using high resolution zoom-in AMR cosmological simulations. In 
the postprocessing, we distinguish between self-shielded gas and unshielded gas. In the 
self-shielded gas, which is optically thick to Lyman continuum radiation, we assume 
pure collisional ionisation for species with an ionisation potential greater than 13.6 eV. 
In the optically thin, unshielded gas these species are also photoionised by the meta- 
galactic radiation. In addition to absorption of radiation from background quasars, we 
compute the absorption line profiles of radiation emitted by the galaxy at the centre of 
the same halo. We predict the strength of the absorption signal for individual galaxies 
without stacking. We find that the Lya absorption profiles produced by the streams are 
consistent with observations of absorption and emission Lya profiles in high redshift 
galaxies. Due to the low metallicities in the streams, and their low covering factors, 
the metal absorption features are weak and difficult to detect. 

Key words: cosmology: theory — galaxies: evolution — galaxies: formation — galax- 
ies: high redshift — intergalactic medium — galaxies: ISM 



1 INTRODUCTION 

Cold gas is thought to flow into massive haloes 10^^ 
at z = 2 — 3 along filaments with velocities of ^ 200kms~^. 
This phenomenon is predicted by simulations and the- 
oretical analysis, where high-z massive galaxies are 
continuously fed by narrow, cold, intense, partly clumpy, 
gaseous streams that penetrate thr ough the shock- 
heate d halo ^as into t he in ner galaxy ('Birnboim & DekeJ 
20031: Keres etHI l2QQ5l : [Pek el & Birnboim 2006s; 

l2008l: iKeres et al.1 l20Q9l: 
Ijohansson, Naab Ostrikerl 
dense, unstable, turbulent 



disc with a bulge and trig g er rapid star formation 
(iDekeL Sari CeverinQ 2009; "Agertz. T^vssier Moord 
20091: LCeverino Dekel k Bournaudl |201q[ 



Krumholz k Burkert boiol; lAgertz. Tevssier k Moord 



201 ll: ICeverino et al.1 [ioT^: ICacciato. Dekel Genelll2012l : 



Sooi; 

Ocvirk, Pichon k Tevssierl 



Dekel et al 
20091 ). 



They 



l2009l : 

form 
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GeneL Dekel k Cacciatol I2OI2I : ICenzel et al.1 120111 ). N- 

body simulations suggest that about half the mass in 
dark-matter halos is built-up smoothly, suggesting that 
the baryons a re also accreted semi-continuously as the 
galaxies grow (iGenel et al.ll2010l ). Indeed, hydrodynamical 
cosmological simulations reveal that the rather smooth gas 
components, including mini-minor mergers with mass ratio 
smaller than 1:10 , brings in about two thirds of the mass 
(jDekel et al.ll2009l ). The m assive, clumpy and star- forming 
discs observed at z ~ 2 dC enzel et al.' '20081: ICenel et al. 
[2OO8; Forster Schreiber et al . 2009, 2011) may have been 
formed primarily via the smooth and steady accretion 
provided by the cold streams, with a smaller contribution 
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by major merger events (lAgertz. Tevssier Moor3 1201 ll : 
ICeverino. Dekel k BQurnaudll2Ql(J )! 



IPettini et alJ (l2002l ): lAdelberger et alJ (l2Q03l ): 
IShaplev et alJ (|2003l ) explored the kinematics of the 
intergalactic medium in the early univer se in ^re a t dept 
using high-resolution spectroscopy. Stei del et al.l (| 19961 ) 
proposed that the predicted inflowing gas could show up as 
redshifted metal absorption lines in the spectra of Lyman- 
break galaxies (LBGs). However, it seems that only a small 
fraction of the LBGs exhibit redshifted absorption features, 
which was interpreted as an indication for the absence of 
cold streams in ha loes with 4 x 10^^ < Mv < lO^^M© 
(jSteidel et aljboiol . hereafter SIO). This has to be reconciled 
with the robust theoretical prediction that the predicted 
filaments are ubiquitous in haloes in the early universe. 
The goal of our paper is to predict the absorption-line 
signatures of the cold flows, for a detailed comparison with 
observations like SIO. 



I Goerdt et al.l (|2Q10l ) used cosmological hydrodynami- 
cal AMR simulations to predict the characteristics of Lya 
emission from the cold gas streams. The Lya luminosity 
in their simulations is powered by the release of gravi- 
tational energy as the gas is flowing with a rather con- 
stant velocity down the potential gradient toward the halo 
centre. The simulated Lya-blobs (LABs) are similar in 
many ways to the observed LABs. Some of the observed 
LABs may thus be regarded as direct detections of the 
cold streams that drove galaxy evoluti on at high z. Ob- 
serv ations s eem to support this pic ture (|Rauch et al ] l201ll : 
Er b. Bogosavlievic Steidel 1201 ih . On the other hand, 
Faucher-Giguere et al.l (l201ol v us ing SPH simu l ations , pre- 



dict lower Lya luminosities than iGoerdt et all (|2010l ). and 
flnd it diflficult to explain the observed LAB luminosities 
with i nflow-driven cooling rad iation alone. This is partly be- 
cause I Faucher-Giguere et al ] (12010) exclude emission from 
very dense gas, and partly because the SPH simulations un- 
derestimate the energy dissipation of the inflowing gas as it 
interacts with itself and with the hot medium. They make 
the point that the gas temperature should also be taken into 
account when distinguishing between the optically thin and 
thick regimes. New AMR simul ations incorporating radia- 
tive transfer seem to confl r m the I Goerdt et al' (I2OIOI) m odel 
(|Rosdahl B laizot 20121 : Kasen et al.l in preparatio"nl ) . In 
particular, Kasen et al.l (|in preparation! ) analvsing tiigh- 
resolution cosmological AMR-hydro simulations with full 
radiative transport for the ionising radiation from back- 
ground and from stars and for the Lya scattering flnd 
that about half the spatially extended Lya luminosity orig- 
inates from fluorescence due to ionising radiation from new 
stars, and the other half from cooling emission due to 
the gravitational energy gain, yieldin g a total lum i nosity 
that is comparable to the flndings of iGoerdt et al.l (|2010l l 
based on crude approximations for the radiative transport 
and without the fluorescence . In a pioneering observation, 
ICantalupo. Lillv k, Haehneltl (|2012l l detected extended Lya 
emission from circum-galactic fllaments of cold gas that 
seem to largely arise from fluorescence excited by the UV 
radiation from an AGN. They detected dark galaxies and 
circum-galactic fllaments that are fluorescently illuminated 
by a quasar. The properties of these Lya sources all indicate 
that the sample is consistent with having much of the Lya 
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Figure 1. Gas density in a simulated galaxy (resolution 70 pc) 
dX z = 2.3 with Mv = 3.5 x IG^^Mo. The colour refers to the 
maximum density along the line of sight. The contours mark n = 
0.1, 0.01 and 0.001 cm~^, respectively. The circle shows the virial 
radius which is 74 kpc. The inflow is dominated by three cold 
narrow streams that are partly clumpy. The density in the streams 
is n = 0.003 — 0.1 cm~^, with the clump cores reaching n ~ 
1 cm~^. The white open boxes tag main features such as streams. 
The same type of points is used to identify these features in figure 



emission originating in fiuorescent reprocessing of quasar ra- 
diation^ 

iBertone k Schavd (|2012l ) investigated the nature and 
detectability of redshifted rest-frame UV line emission from 
the IGM at 2 < 2; < 5. They used simulations to create 
maps of a set of strong rest-frame UV emission lines. They 
conclude that several species (e.g., C III, C IV, Si III, Si IV, 
and O VI) of those from the high-redshift IGM will become 
detectable in the near future. They also claim that lower 
ionisation lines provide us with tools to image cold accretion 
fiow s as well as cold outflowin g clouds. 

IVaZ de Voort et al 1 (l2012l ). using SPH simulations, flnd 
that nearly all of the Hi absorption arises in gas that has 
remained fairly cold, at least while it was extragalactic. In 
addition, the majority of the H i is rapidly falling towards a 
nearby galaxy, with non-negligible contributions from out- 
flowing and static gas. They identify carefully the environ- 
ment of the Hi absorbers and conclude with this informa- 
tion that cold accretion flows are critical for the success of 
simulations in reproducing the observed rate of incidence of 
damped Lya and particularly that of Lyman limit systems. 
They propose that cold accretion flows exist and have al- 
ready been detected in the form of high column density H I 
absorbers. 
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Figure 2. Distribution by volume of the neutral hydrogen fractions xm and total hydrogen gas densities nu in the circumgalactic 
environment of our simulated galaxy (see text). The left panel shows the distributions for the full radiative transfer calculation including 
the metagalactic UVX field plus internal sources presented by Fll. The middle panel is for the radiative transfer calculation, but with 
UVX only. The right panel are results employing a pure density criterion of tir, shield = 0.01 cm~^ for the boundary between self-shielded 
and optically thin cells. The units of the colour-bar axis are unity per unit area in the (logio -^Hl, logio ^h) plane. Contour lines are at 
0.0001, 0.001 and 0.01 respectively. Our simplified ionisation model provides a useful crude approximation for the UVX radiative transfer 
results (see text). 



iKacprzak et al. 1 (l2010>) analysed synthetic Mgll absorp- 
tion spectra from simulated galaxies ai z — 1 and they ten- 
tatively concluded that most of the Mgii absorption arises 
in filaments and tidal streams of the galactic halo. This ma- 
terial is infalling with velocities around 200 km s~^. They 
claim that large equivalent- width (EW) Mg ii are associated 
with outflows. 



iFaucher-Giguere k Kered (l201lh. iKimm et al.1 (l201lh . 
IStewart et al. ~ (|201lh and iFumagalli et al .1 (I2OIIL hereafter 



Fll) estimated from simulations the covering factor of the 
cold streams at ^ = 2 — 4 for absorption from a back- 
ground sou rce. They found that it is small and decreasing 
with time. IStewart et al ] (l201lh and Fll argue that as a 
galaxy switches from cold- to hot-mode accretion with de- 
creasing redshift the reduced cold gas accretion naturally 
results in a suppression of cool circumgalactic gas. This 
transition occurs when the halo is mass ive enough to sup- 
port a stable shock near the virial radius (iBirnboim k Dekell 
I2OO3I : iKeres et al.ll2005l : IPekel k Birnboimll2006l ). Fll dis- 
cover that within ^ 500 Myr of reaching this threshold mass, 
the covering fraction drops. Because the cold-stream cover- 
ing factor is small compared to the order unity covering 
fraction expected for galactic winds, the cold streams are 
generally overwhelmed by outflows in absorption spectra. 
However, Fll show that there is enough cross section in the 
cold flows that could make them detectable in Lya absorp- 
tion in surveys of background sources. Note the differences 
with the observations of galaxies at zero impact parame- 
ter, i.e. observations in the transve rse direction, using e.g. 
quasars. IVa^ de Voort et al ] (I2OI2I ) and Fll show that in 
fact there is a significant contribution to the population of 
absorption line systems that comes from the streams. That 
is, outflows are prominent, but inflows should be visible as 
well. 

iKimm et al.l (|201lh and Fll investigated whether the 
cold gas streams are expected to be detectable in low- 
ionisation metal absorption lines, such as C 11 (1334 A) or 



Sill (1260 A). They find that the metal absorption signal 
by the interstellar medium of the central galaxy itself is so 
deep and so broad in velocity space that it swamps the metal 
signal from the filamentary gas, which is rather metal poor. 
A cold filament might be detectable in metal lines if it lies 
precisely along the line-of-sight, but this would be rare. 

IShen et al.l (|2012l ) used SPH simulations adopting a 
blastwave scheme for supernova feedback which, in com- 
bination with a high gas density threshold for star forma- 
tion, has been shown to generate large-scale galactic out- 
flows. Their strong outflows avoid the incoming dense nar- 
row streams, and flnd their way out in a wide solid angle 
through the dilute medium between the streams. Showing 
that the cold inflows are not affected by the outflows. They 
generate synthetic spectra by drawing sightlines through the 
simulated CGM at different galactocentric impact parame- 
ters, and compared the theoretical interstellar absorption 
line strengths with the observations. They flnd that they 
are in broad agreement with those observed at high-redshift 
by SIO. 

As in the observations of SIO, we consider separately 
absorption features of radiation emitted from background 
sources and from central sources, namely along lines of sight 
with flnite and with vanishing impact parameter with re- 
spect to the galaxy at the centre of the halo that hosts 
the absorbing gas. A theoretical study of absorption in the 
central-source geometry has not been attempted so far. The 
same is true for a statistical study of absorption signal in in- 
dividual galaxies without stacking. We focus on predicting 
absorption line proflles that correspond to the measurements 
of SIO. We refer to the same ten different lines and mimic 
their observational procedure by averaging data for several 
simulated galaxies and many lines of sight. As in Fll, we em- 
ploy high-resolution AMR simulatio ns that are described in 
ICeverino. Dekel k Bourn^ l\201(t ). The spatial resolution 
in these simulations is 20 times better than in the Horizon 
MareNostrum simulation (|Ocvirk. Pichon k Tevs"sie3l2008l ) 
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Figure 3. Hammer equal area projections of Hi column densities in our fiducial simulated galaxy (resolution 70 pc, z = 2.3, Mv = 
3.5 X IQ-^-^Mq). The upper panel is integrated from 0.3 to 1.0 J?v, the lower panel from 0.1 to 0.3 Rv The contour lines indicate 10^° 
and lO-'^^cm"^. The narrow major streams with column densities up to ~ 10^^ cm~^ can be seen, as well as the vast majority of the 
stream-free sky (as seen from the central galaxy) having only column densities of < 10-*^^ cm"-^. The white open boxes in the upper panel 
identify the main features such as the streams. The same type of points is used to identify these features in figure [T] The reader should 
be aware of the fact that these plots are theoretical and not observed in nature since one always looks down to the very centre at r = 0. 



used bv lKimm et al.l (| 20 111 , they also used higher resolution 
simulations for verification purposes). Unlike in Fll, we do 
not employ a full radiative transport analysis, and instead 
account for self-shielding by a simple density criterion that 
approximates the radiative-transport results. The compu- 
tations include explicitly highly ionised species (e.g. C iv) 
created in the warm highly ionised intergalactic medium (as 
opposed to C II and others that are created in the cold low- 
ionised neutral medium). 

Our paper is organised as follows. In section [2] we in- 
troduce the simulations used for this analysis. In section [3] 
we explain the computation of the absorption features. In 
section [4] we show the results for the central geometry. In 



section [5] we show the results for the background geometry. 
Finally, in section [6] we draw our conclusions. 



2 THE SIMULATIONS 

We use three simulated galaxies from a suite of simulations, 
employing Eulerian adaptive mesh refinement (AMR) hy- 
drodynamics in a cosmological setting. These are zoom-in 
simulations in dark-matter haloes with masses ~ 5x 10^^ 
at z = 2.3, with a maximum resolution of 35 — 70 pc in phys- 
ical coordinates ([Ceverino. Dekel BournaudI |201Q| . here- 
after CDB). Figure [T] shows a density map of one of the 
three galaxies from the simulation, which serves as our fidu- 
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cial galaxy. It demonstrates the domi nance of typically thre e 
co-planar narrow cold streams (see IPanovich et al.l l20l3 v 
originating outside the virial radius along the dark-matter 
filaments of the cosmic web, penetrating into the discs at 
the halo centres. The streams consist of a smooth compo- 
nent and clumps with a spectrum of sizes. The typical den- 
sities in the streams are in the range n = 0.003 — 0.1 cm~^, 
and they reach n = 1 cm~^ near the central disk and at the 
clump centres. Some of those clumps are satellites with dark 
matter haloes. 

The CDB simula tions were run with the ART ( adap- 
tive refinement tree; iKravtsovT Klvpin Khokhlovl Il997l : 
iKravtsovl l2003l ) code. It incorporates relevant physical 
processes for galaxy formation, including gas cooling, 
photo- ionisation, heati ng, star formation, metal e nrichment 
and stellar feedback (|Ceverino Klvpinl l2009l ) . Cooling 
rates were computed for the given gas density, tempera- 
ture, metallic i ty an d UV background (based on cloudy, 
iFerland et al.l Il998h . Cooling is assumed at the centre 
of a cloud of uniform d ensity with thickness 1 kp c 
(ICeverino-Rodriguej l2008l : [Robertson Kravtsovl l2008l ). 
Metallicity dependent, metal- line cooling is included, assum- 
ing a relative abundance of elements equal to the solar com- 
position. The code implements a "constant" feedback model, 
in which the combined energy from stellar winds and super- 
nova explosions is released as a constant heating rate over 
40 Myr. This is the typical age of the lightest star that ex- 
plodes as a type-II supernova. Photo-heating is also taken 
into account self-consistently with radiative cooling. A uni - 
form UV background based on the lHaardt Madaul (Il996l ) 
model is assumed. Local sources are ignored. 

This code has a unique feature for the purpose of simu- 
lating the detailed structure of the streams and the gravita- 
tional instability in the disk. It allows gas cooling to well be- 
low lO^K. This enables high densities in pressure equilibrium 
with the hotter and more dilute medium. A non- thermal 
pressure floor has been implemented to ensure that the 
Jeans length is resolved by at least seven resolution elements 
preven ting artificial f ragmentation on the smalles t grid 
scale dT ruelove et al." 1997; Rob ertson Kravtsovl l2008l : 
ICeverino , Dekel & Bournaud 2010). The pressure floor is ef- 
fective in the dense {n > 10 cm" and cold (T < lO^i^) 
regions inside galactic disks. Most of the absorbing gas in 
question, which is far outside the disk, is not affected by 
this pressure floor. 

The equation of state remains unchanged at all densi- 
ties. Stars form in cells where the gas temperature is below 
lO^K and its density is above the threshold n = 1 cm""^ 
according to a stochastic m odel that is roughly consistent 
with the iKennicuttI (|l998l ) law. The ISM is enriched by 
metals from supernovae type II and type la. Metals are 
released from each stellar particle by SNII at a c o nstant 
rate for 40 Myr since its birth. A iMiller Scalol (|l979l ) 
IMF is assumed. This p rocedure matches the results of 
IWooslev Weaverl (|l995l ). The metal ejection by SNIa as- 
sumes an exponentially declining SNIa rate from a max- 
imum at 1 Gyr. The code treats the advection of metals 
self-consistently and it distinguishes between SNII and SNIa 
ejecta (Ceverino- Rodriguez 2008). 

The dark matter particle mass is 5.5 x lO^Mo.The mmi- 
mum star particle mass is lO^M©. The initial conditions for 
the CDB simulations were created using a low- resolution 
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Figure 4. Cumulative Lya sky covering fractions fc higher than 
column density A^hi- It is integrated from 0.3 to 1.0 i?v for all 
gas inflowing and outflowing. For the blue dashed line we use 
our three simulated galaxies (resolution 70 pc) at z = 2.3 with 
Mv = 3.5 X 10^1 M0 and i^v = 74 kpc, for the red solid hue we 
use two simulated galaxies at z = 1.38 with Mv = 1.38 x 10^^ Mq 
and Rv = 150 kpc. The observational visibility limit at 10^^ cm"^ 
is indicated by the vertical line. We see very low (< 0.3) covering 
fraction over almost the whole range of column densities and a 
a sharp jump of the covering fraction at the visibility limit. The 
higher redshift galaxy has a considerably higher covering fraction 
over the whole range of column densities. 

cosmological A/'-body simulation in a co-moving box of side 
29 Mpc. Its cosmological par ameters were motivated by 
WMAP5 (iKomatsu et al.ll2009l ). The values are: 0^ = 0.27, 
0.73, Qh 0.045, h 0.7 and as 0.82. At z = 1, 
three haloes of Mv — 10^^ Mq each have been selected, 
avoiding haloes that were subject to a major merger at that 
time. The three halo masses at z = 2.3 are 3.5, 4, 6x 10^^ M©, 
their virial radii are around 74 kpc and they end up as 
(3 — 4) x 10^^ M0 haloes today. Around each halo, a con- 
centric sphere of radius twice the virial radius was selected 
for resimulation with high resolution. Gas was added to the 
box following the dark matter distribution with a fraction 
/b = 0.15. The whole box was then re-simulated, with re- 
fined resolution only in the selected volume about the re- 
spective galaxy. 



3 COMPUTING THE IONISATION STATES 

A reliable estimate of the ionisation state of the gas is es- 
sential for any study of absorption line systems. Ideally, the 
determination of the ionisation states of the gas should be 
coupled to the hydrodynamic calculations, but due to the 
numerical complexity of the problem this type of calcula- 
tion remains computationally expensive for high resolution 
simulations especially at lower redshift s. In order to com- 
pute them via post-processing we take the densities, tem- 
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Figure 5. Same cumulative sky covering fractions as in figure @]this time for metal lines. In the right panel we use our three simulated 
galaxies (resolution 70 pc, z = 2.3, Mv = 3.5 x 10^^ Mq and i?v = 74kpc), in the left panel we use two simulated galaxies at z = 1.38 with 
Mv = 1.38 X 10^^ Mq and = 150 kpc. The observational visibility limit at lO^^cm"^ is indicated by the vertical line. Together with 
figure 0] this plot shows that a full sky map of column densities exhibits a few very high peaks corresponding to the streams extending 
only over a fairly tiny solid angle and a majority of very low column density solid angle corresponding to the space without streams. 
This effect is stronger for the low redshift panels. 



peratures and metallicities from the simulation. We assume 
a primordial helium mass fraction Y — 0.24, corresponding 
to a helium particle abundance of 1/12 rela tive to hydrogen . 
For the heavy elements we assume the iAsplund et alJ (|2Q09l ) 
solar photosphere pattern. 

We determine whether a given simulation cell is "self- 
shielded" , i.e. optically thick to Lyman continuum radiation, 
via a simple density criterion. Cells with total hydrogen ex- 
ceeding rishieid = 0.01 cm~^ are assumed to be self-shielded. 
For lower densities the cells are assumed to be optically thin. 
We find that setting nghieid = 0.01 cm~^ gives the best agree- 
ment with a radiative transfer calculation as we demonstrate 
in figure [2] (see discussion below). For a furt her discuss ion 
of the density criterion, see the discussion bv lCoxl (|2005l l or 
the appendix A2 of Fll. W e calculate the atom ic and ionic 
fractions XAi using cloudy (|Ferland et al.|[l998l ), where XAi 
is defined as the fractions of element A in ionisation state 
i. The specific species we are interested in are Hi, Oi, Cii, 
C IV, Si II and Si iv. These are the most important ions that 
are visible and detected in optical spectra at z = 2 — 3 and 
discussed in great detail by SIO. We also consider Mgii and 
Fell since they are assumed to be very suitable for future 
observational programmes. Table [T] lists the wavelengths and 
oscillator strengths of the various absorption lines we con- 
sider. 

In our CLOUDY computations we distinguish between 
self-shielded and optically thin cells as follows. For optically 
thin cells, we ass u me t hat the gas is exposed to the full 
iHaardt k Madaul (|l996h UV to X-ray (UVX) metagalac- 



ion 


Ao [A] 


fx 


7A [s-i] 


Lyo; 


1215.6701 


0.4164 


6.265 X 10^ 


Cii 


1334.5323 


0.1278 


2.870 X 10^ 


Sill 


1260.4221 


1.115 


2.533 X 10^ 


Oi 


1302.1685 


0.04887 


5.750 X 10^ 


Sill 


1526.7066 


0.1155 


1.960 X 10^ 


Sill 


1304.3702 


0.09345 


1.720 X 10^ 


Civ 


1548.195 


0.1908 


2.654 X 10^ 


Si IV 


1393.755 


0.514 


8.825 X 10^ 


Civ 


1550.770 


0.09522 


2.641 X 10^ 


Si IV 


1402.770 


0.2553 


8.656 X 10^ 


Mg 11 


2796.352 


0.6123 


2.612 X 10^ 


Fell 


2382.765 


0.3006 


3.100 X 10^ 



Table 1. Basic parameters of the various absorption lines, we 
quote wavelength Aq, oscillator strength fx and the damping 
width '^x which is the su m over t he spontaneous emission coeffi- 
cients. Data is taken from lMortonI (^991), exce pt for the /a - values 
for the three Sill lines, which are from lDufton et al.l (|l983[ ). 



tic background field at the appropriate redshift. We do not 
include any radiation from "local sources". The ionisation 
state is then controlled by the combined effects of electron- 
impact collisional ionisation and photoionisation for the full 
spectral range of the background field. For shielded cells, we 
truncate the radiation field and include only photons with 
energies below the Lyman limit. Hydrogen and other species 
with ionisation potentials (IP) greater than 13.6 eV are then 
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Figure 6. Example Lya absorption line profiles seen from a single direction in the central source geometry integrated from 0.3 to 1.0 
Rv. Positive velocities are infiowing into the galaxy and negative velocities are out of the galaxy. The letters correspond to the indicated 
position of the Hammer projection from the upper left panel of figure [S] Values for neutral hydrogen column density and EW are quoted 
in each panel. Uppermost panel: various modes of infiowing streams, intermediate upper panel: different examples of outfiowing material, 
intermediate lower panel: low absorption profiles and lowermost panel: extraordinary cases which are all showing clear signatures of 
infiows. Some of the lines are saturated like A, C, K and L whereas others are not. 
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Figure 7. Same figure as figure (6] panel D (the panel with the deepest signal) this time for selected metal lines providing also the 
deepest signal, namely: Cll (1334 A), Sill (1260 A), Mgll (2796 A) and Fell (2383 A). Panel D is fully saturated in Lyo; but still shows 
peak absorp t ion li ne depths of 0.2 and EW of ~ 0.1 A for these metal lines. The C ll panel of this plot disagrees with the results of 
iKimm et al.l (|201lL their figure 2, lower panel) due to different line computing algorithms. 



produced by collisional ionisation onljQ. For the shielded 
cells we neglect photoionisation by penetrating X-rays. So 
for pure collisional ionisation in s hielded cells the hydroge n 
becomes ionised for T > 2 x 10"^ K (jCnat SternberdlioOTi ). 
However, "low-ions" with IPs below 13.6 eV continue to be 
produced by photoionisation as well. Thus, for a given cell, 
the ionisation state is determined by its total hydrogen den- 
sity flu and temperature T, and the adop ted radiation field, 
which is either the "full" or "truncated" iHaardt Madaul 
(| 19961 ) field, depending on whether the cell is characterised 
as shielded or not. 

To verify our procedure for distinguishing between 
shielded and optically thin cells, we compare our results for 
the hydrogen ionisation states with the more detailed ra- 
diative transfer computations presented by Fll. In figure [2] 
we plot the distributions by volume of the neutral hydro- 
gen fractions xhi and total hydrogen gas densities nu in 
the simulation cells, for three different computations of the 
ionisation structure, all for the same simulated galaxy. The 
left-hand panel shows the neutral fraction versus density dis- 
tributions for the detailed radiative transfer computation of 
Fll. This includes the combined effects of the external UVX 
and internal sources of photoionising radiation. The middle 
panel shows the Fll results for UVX only. Cells with den- 
sities higher than 0.1 cm~^, which usually are close to the 
disc or the satellites, are not included in figure [2l The right- 
hand panel shows our results for xm versus nn using our 
simple criterion of assuming the cells are fully shielded for 
> '/T'H, shield = 0.01 cm~^. Again, in the post-processing 
we assume that the hydrogen is collisionally ionised in the 
shielded cells, and collisionally plus photoionised in the op- 
tically thin cells. It is apparent that this simplified model 
provides a useful crude approximation for the results of the 
UVX radiative transfer (middle panel). The only significant 
disparity occurs near the transition density of 0.01 cm~^ and 
neutral fraction of xhi ~ 0.1, where the ionisation states 
appear to be sensitive to the details of the radiative trans- 
fer computations. At the end of section [5] we will compare 
our "final product" (figures [24] and [25)) to similar results al- 



^ The reader should note that all our species listed in Table [T] 
have IPs greater than 13.6 eV. 



ready published in the literature computed with full radia- 
tive transfer including ionising radiation from local sources 
(figure 13 of Fll). The differences are small, justifying our 
simplifying model. 



4 CENTRAL SOURCE 

In this section we consider the Lya and metal line absorption 
that occurs as UV light emitted by the central galaxy is ab- 
sorbed by gas in the circum-galactic environment. As defined 
by SIO, the circumgalactic medium is situated in the spher- 
ical zone from just outside the galactic disk to around the 
virial radius R^. Observations of absorption against the cen- 
tral galaxy itself have the advantage of being able to discrim- 
inate between inflows and outflows because the absorptions 
may be assumed to occur in foreground material only. Radi- 
ation emitted or scattered from behind the galaxy is blocked 
by the galaxy itself. However, such absorptions do not pro- 
vide spatial information about the distance from the galaxy 
centre as they are all by definition at an impact parameter 
b — from the galaxy centre. SIO employed this technique 
of observing absorptions against the central galaxy by stack- 
ing a sample of 89 galaxies with z — 2.3 di 0.3 using both 
rest-frame far-UV and Ha spectra, to investigate the kine- 
matics of the gas flows in the circumgalactic regions. Here 
we predict the spectral line profiles and equivalent widths 
(=EW) produced in absorption against the central galaxy, 
for a direct comparison to the SIO observations (as presented 
in their section 4). 

The absorption-line profiles depend on the column den- 
sities along the line-of-sight. We compute the column den- 
sities along radial rays from an innermost radius ri (usually 
around the galactic disk) to the virial radius Rv The col- 
umn A^Ai(0, 0) of element A in ionisation state i, at angular 
position (0, 0) is calculated by 



A^Ai(0,^) 



XAi(r) nA{r) dr, 



(1) 



where ua (f) is the total gas density of element A at position 
r = (0, 0, r) and xai is the ionisation fracti on. In figure [3] we 
show two Hammer equal area projections (ISnvderlll993h of 



© 2012 RAS, MNRAS OOO.fTti26l 



Absorption in cold streams feeding high-z galaxies 9 




0.0001 0.001 0.01 . 0.1 1 -300 -200 -100 100 200 300 

Wo [A] outflow Vr[km/s] inflow 



Figure 8. Hammer equal area projections of EW (left panels) and of EW-weighted velocity (right panels) in our fiducial galaxy. The 
upper four panels show Lyo; whereas the lowermost panels show C ll. In the uppermost as well as in the lowermost panels the data were 
computed between 0.3 and 1.0 R^, whereas in the middle panel the data were computed between 0.1 and 0.3 Rv Positive velocities 
are infiowing into the galaxy (shown in red) and negative velocities are out of the galaxy (shown in blue). The contour lines indicate 
EWs of 1.0, 0.1, 0.01, 0.001 and 0.0001 A or velocities of +200, and -200 kms-\ respectively. In the upper left panel the positions 
of the example line profiles from figure [6] are indicated by the circles having letters attached. Comparing right and left panels gives a 
comprehensive picture of infiows and outfiows. The velocity pattern of Lyo; and Cll is the same even if their EWs are very different. In 
the panels spanning the outer radius range the narrow major streams are clearly seen as well as the vast majority of the "stream- free" sky 
(as seen from the central galaxy). In the middle panel representing a narrow radius range very close to the central galaxy this bi-modality 
is not as obvious as in the other two cases. The radial velocities show different behaviour in the two different radius regimes, the Lya 
and C II panels of the same radius regime on the other hand are very similar. 



Hi column densities for our fiducial galaxy (z = 2.3, Mv = 
3.5 X 10^^ Mo, i?v = 74kpc). In the upper panel the column 
densities are integrated from 0.3 to 1.0 whereas in the 
lower panel they are integrated from 0.1 to 0.3 R^. Excluding 
the very inner sphere is somewhat artificial since in practise 
absorption will occur from close to the very centre at r = 0. 
However, these plots are useful in giving an impression of the 
distribution of column densities along different sight-lines. In 
the upper panel the main stream features are identified and 
tagged to the same features shown in figure [1] by white open 
symbols. The same symbols types in figures [T] and [3] indicate 



identical features in the galaxy but for different geometries. 
In figure [3] the narrow major streams with column densities 
up to ^ 10^^ cm~^ can be seen, as well as the vast majority 
of the "stream-free" sky (as seen from the central galaxy) 
having only column densities of < 10^^ cm~^. 

To quantify the column density distribution, we com- 
pute the cumulative fraction fc of the sky (as viewed from 
the central galaxy) that is covered with column densities 
greater than A^Ai- The covering fractions for hydrogen and 
the various metal ions are plotted in figures |4] and \5\ These 
plots include inflowing and outflowing gas. The data in the 
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blue dashed line of figure |4] as well as in the right panel of fig- 
ure [5] is averaged over our three simulated galaxies at z = 2.3 
with Mv = 3.5 X lO^^M©. For better comparability with fu- 
ture observations we also use a second set of data consisting 
of two galaxies at 2; = 1.38 with Mv = 1.4 x 10^^ M© (red, 
solid line). The plot shows that the whole sky, as seen from 
the central galaxy, has a minimum A^hi of 10^^ cm~^ and 
above this value fc decreases strongly from 1.0 to 0.3 at 10^^ 
cm~^. The Hi covering fraction then decreases more slowly 
down to ^ 1% at 10^"*^ cm~^. The fc for the metals are of 
course much lower. The highest metal columns are usually 
around 10^^ cm~^ and at 10^ cm~^ the covering fractions 
vary between 0.10 (Fell) and 0.75 (Civ). This verifies our 
initial qualitative impression that a full sky map of column 
densities exhibits a few very high peaks corresponding to the 
streams extending only over a fairly tiny solid angle and a 
majority of very low column density solid angle correspond- 
ing to the space without streams. This effect is even more 
impressive for the low redshift panels. Also there is a de- 
crease of the covering fraction with decreasing redshift. 

Given the computed column densities we can construct 
simulated absorption line profiles. We start with a single 
line of sight (i.e. a single pixel in figure [3|). Along a given 
sight-line the gas has a varying density, temperature and ra- 
dial velocity as function of radial position r from the central 
galaxy. The convolutions of the different densities and radial 
velocities in the gas along the line of sight are the main in- 
gredients to compute an absorption line profile. This is done 
as follows: the radial velocity offset Aw relative to central 
source at rest of all the gas is measured. We assume Voigt 
profiles with a thermal Doppler broadening parameter 



2k T 

ruA 



(2) 



where k is the Boltzmann constant, T is the temperature 
of the gas and mA is the mass of the element A. So for 
angular position (0, ^) we can compute the optical depth 
Tu {(i>,0, Aw) at the vel ocity offset Aw as (cf. equation 10.24 
of lBohm-Vitenselll990l ) 



T^[(^,0,Aw) 



fx A, 



H 



7a Ao Aw ■ 



b{f) 
- v{r) 



47V b(r) 



b{r) 



dr. 



(3) 



where e is the electron charge, me is the electron mass, c is 
the speed of light, Ao is the transition wavelength, nA{r) is 
the gas density of element A at position f= (0, 0, r), XAi is 
the ionisation fraction of element A in state v is the veloc- 
ity of the gas, fx is the oscillator strength of the absorption 
line and 7a is the sum over the spontaneous emission co- 
efficients or the damping width. We took these fx and 7a 
values from iMortonI (|l99lh . The Voigt profile is 



H{a,u) = - 



exp(-^ J 
(u - yY + a" 



(4) 



where a is the ratio of the damping width to the Doppler 
width and u is the offset from line centre in units of Doppler 
widths. Knowing Tv{Aw) it is now possible to compute an 
absorption line profile I (Aw) for a given direction which is 
simply the function 
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Figure 9. Cumulative EW sky covering fractions fc higher than 
Wo for a central source integrated from 0.3 to 1.0 Rv for Lya. 
Please note the different x-axis scaling with respect to figure IIUI 
In solid red we show the average of our three simulated galaxies 
(resolution 70pc, z = 2.3, Mv = 3.5 x IO^^Mq, = 74kpc). 
Whereas in dashed blue we show the average of two simulated 
galaxies at z = 1.38 (Mv = 1.38 x R^ = 150 kpc). 

At higher redshift there is a higher covering fraction of high EW 
visible, implying that there are more pronounced streams present. 



I (Aw) = exp[— Tjy(Ai(;) 



(5) 



For a fair comparison to the observations done by SIO, we 
mimic a Gaussian point spread function. This will be done 
from now on up to the end of this section. It is in contrast to 
figures [3l|4] or [5] where no point spread function was applied. 
It has a beam-size (= FWHM 77 of the Gaussian) of 4 kpc. It 
is done by splitting up a cylinder with a radius of three times 
the beam-size into as many parallel fibres as the resolution 
permits, determining the absorption line profile for every 
individual fibre and then computing a Gaussian weighted 
average absorption line profile from all fibres. Also from now 
on up to the end of this section (except for figures [TT] and [12]) 
we degrade our velocity space resolution down to 50 km/s, 
to match the velocity resolution used in the corresponding 
figures of SIO. 

A selection of the large variety of different possible ab- 
sorption line profiles depending on the respective viewing 
angles are shown in figures [6] (Lya) and [7] (selected metal 
lines). One can see various modes of infiowing streams, dif- 
ferent examples of outfiowing material as well as very shal- 
low absorption line profiles corresponding to a line of sight 
with very low column density. In Lya some of the lines are 
saturated like A, C, K and L, whereas the others are not. 
Panel D which shows the strongest lines is fully saturated in 
Lya but still shows peak absorption line depths of 0.2 and 
EW of ^ 0.1 A for selected metal lines. The reader should 
note that some of the weaker lines (panels G or H for ex- 
ample) might be completely erased by noise and the ISM 
component at v = 0. 
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Figure 10. Same cumulative EW sky covering fractions fc as in figure [9] this time for the metal lines. Note the different x-axis scaling 
with respect to figure O In the upper panels we show our three simulated galaxies (resolution 70 pc, z = 2.3, Mv = 3.5 x 10^^ Mq, 
Rv = 74kpc). Whereas in the lower panels we use two simulated galaxies at z = 1.38 with Mv = 1.38 x IO-'^^Mq and = 150 kpc. 
The CiV and Si IV lines in the middle panels clearly stand out. Those lines can be produced in hot medium whereas all the others can 
only be produced in cold gas only. So C IV and Si IV reach covering fractions as high as unity. The low ionised metal lines (left and right 
panel) have higher peak values of EW by up to one order of magnitude. 



The leftmost p anel of fi^ure[7 is com parable to the lower 
panel of figure 3 in lKimm et alj ( 201lh . They get for Cii a 
maximum line depth (MLD) of ^ 0.9 with a FWHM ry of 
650 km s~^. The differences stem from the fact that they 
used a lower res olution simulation, the Horizon MareNos- 
trum simulation (|Ocvirk, Pichon Tevssie'3l2008l V Further 
on they used higher mass haloes (Mv > 10^^ M©) and a 
very different prescription for the Gaussian velocity distri- 
bution. They obtain their dispersion from the neighbouring 
26 cells instead of the algorithm from equations ([2|) to Q. 
Some statistical properties of the ensemble of all possible 
example line profiles will be shown later. 

Each of these absorption line profiles has a certain to- 
tal absorption which can be expressed as a single number, 
the "rest equivalent width" (EW). This is usually given in 
wavelength units and it is defined as: 



Wo = 



Ao 



{1 — exp[— Tiy(i(;)]} dw 



-c 



[1 - I{w)] dw. 



(6) 



An EW can now be calculated for every angular position 
on the sky. With these numbers we can plot a whole full 
sky map of EW. This is done in the left panels of figure [51 
There we show the resulting Hammer projections of EW for 
our fiducial galaxy. The upper two panels show Lya whereas 
the lowermost panel shows Cii. In the uppermost and low- 
ermost panels the EWs were computed for the gas in the 
radius range between 0.3 and 1.0 i?v, whereas in the mid- 
dle panels the EWs were determined for gas in a somewhat 
narrower and innermore radius range between 0.1 and 0.3 
i?v. In the upper left panel the positions of the example 
line profiles from figure \6\ are indicated by the circles with 
the letters attached to them. In the panels spanning the 
outer radius range one can again clearly see the narrow ma- 
jor streams having EWs of > lA (Lya) or > 0.01 A (Cii). 
One can also see the vast majority of the "stream-free" sky 
(as seen from the central galaxy) with EWs of ^ 0.01 A 
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Figure 11. Fraction h^ip of all possible example line profiles in 
the central geometry whose EW is higher than Wq and whose line 
centre is at a velocity offset indicating an inflow at least as fast 
as Aw. Here for the Lyo; line. The example line proflles are inte- 
grated from 1.0 i?v down to 0.3 Rv Positive velocities are inflow- 
ing into the galaxy and negative velocities are out of the galaxy. 
Contour lines are at 0.1, 0.01, 0.001 and 0.0001 respectively. Note 
the different y-axis scalings with respect to flgure [T2] This flgure 
shows that the probability of detecting an inflow which is flowing 
in with at least 150 km s"-*^ with a signal of at least 3.9 A in a 
single observation of a single galaxy without averaging is ~ 4%. 



(Lya) or < 10~^ A (Cli). In the middle panel represent- 
ing a narrow radius range very close to the central galaxy 
(0.3i?v > r > 0.1i?v) this bi- modality is by far not as obvi- 
ous as in the other two cases. This is due to the fact that at 
this distance from the galaxy centre the inflowing streams 
slowly dissipate into a region with more complex gas motions 
("messy region"). Since the EWs are significantly lower for 
the C II line and even lower for the other lines, as we will 
see later, we do not show Hammer projections of EW of 
other metal lines than C ii. We do show Hammer projections 
of EW- weighted velocity in the right panels of this figure. 
Positive velocities are inflowing into the galaxy (shown in 
red) and negative velocities are out of the galaxy (shown in 
blue). The contour lines indicate velocities of +200, and 
-200 kms"^ Re gions with very high velocity outflows have 
a low EW, while velocities close to systemic velocity are at 
very high EW. This agrees very well with the high EW seen 
by SIO at Vr = 0kms~^. The velocity pattern of Lya and 
Cii is the same even if their EWs are so different. Compar- 
ing the upper two panels one can see that the distribution 
of velocities in the two different radius regimes varies. The 
term "outflow" can be misleading here. In our simulations 
some of the material that is receding from the centre of the 
galaxy is not necessarily diffuse gas driven by some kind 
of feedback processes. It could also be a clump or satellite 
galaxy which is merging with the central galaxy. It just had 
its first flyby with the central object and is therefore now 
"outflowing" . 

To present the statistics from these maps we determine 
the sky covering fractions. In figures [9l and [TOl we show the 
fraction fc of the sky as seen from the central galaxy that is 



covered with absorption line profiles which have EW higher 
than Wq. Here we integrate in the central source geome- 
try from 0.3 to 1.0 R^. In figure [9] this is done for Lya 
whereas in figure [10] this is done for the metal lines (note 
the different x-axis scaling in both figures). We use our 
three simulated galaxies (resolution 70 pc) at 2; = 2.3 with 
Mv = 3.5 X 10^^ M© and R^ = 74kpc. These are shown in 
solid red in figure [9] and on the upper panels of figure [101 
Again for better comparability with future observations we 
also use a second set of data consisting of two galaxies at 
z = 1.38 with Mv = 1.4 x lO^^M© and R^ = 150 kpc. These 
are shown in dashed blue in figure [9] and on the upper pan- 
els of figure [TOl We concentrate on the lines studied in great 
detail in SIO which are: Lya (1216 A), C 11 (1334 A), Ol 
(1302 A), Sin (1260 A, 1304 A and 1526 A), Civ (1548 A 
and 1550 A) as well as Si IV (1393 A and 1402 A). In ad- 
dition to these ten lines observed by SIO we also use here 
seven new ones: Mgii (2803 A and 2796 A) and Fell (2600 
A, 2383 A, 2344 A, 2587 A and 2374 A). These lines have 
not been observed in this geometry or redshift range yet, but 
may be very suitable for future observational programmes 
at 2; 1. For Lya one can see that at higher redshift there 
is a higher covering fraction of large EW, implying that the 
streams are much more pronounced. The C IV and Si IV lines 
in the middle panels of figure [TOl clear Iv stand out since they 
reach covering fractions as high as unity but for low column 
densities only. This is because these lines can be produced in 
the hot and diffuse medium whereas all the others can only 
be produced in cold gas only. Another interesting feature for 
the metal lines is that low ionised metal lines have higher 
peak values of EW by up to one order of magnitude, but low 
overall covering fraction. This is because low-ionised ions are 
more abundant in high-density, cold regions with low cover- 
ing fraction, like streams. 

As already mentioned earlier we want to show statistics 
of how many absorption lines show signatures of inflow ver- 
sus how many absorption lines show signatures of outflow 
and how strong are those. This is done in figures [TTI and [T2l 
These figures are intended to give the likelihood of a detec- 
tion of a cold stream while looking at a single galaxy from a 
single direction without averaging. They show the fraction 
help of all possible example line profiles in the central geom- 
etry whose EW is higher than Wo and whose line centre is at 
a velocity offset indicating an inflow at least as fast as Aw. 
The example line profiles are integrated from 1.0 i?v down 
to 0.3 i?v Positive velocities are inflowing into the galaxy 
and negative velocities are out of the galaxy. In these two 
figures no velocity degrading is used. The Lya panel shows 
that the probability of detecting an inflow which is flowing 
in with at least 150kms~^ with a signal of at least 3.9 A 
in a single observation of a single galaxy without averaging 
is ^ 4%. In C II one sees an inflow > 150kms~^ with an 
EW > 0.17 A in 0.4% of all observations. For Mgii one sees 
an inflow > 150 km s"^ with an EW > 0.2 A in 1.3% of all 
observations. These values should be achievable by future 
observations. Since it is the likely scenario that the inflow- 
ing streams and the wide-angle outflowing gas avoid each 
other to a great extent the streams as measured in our sim- 
ulations where the outflows are a factor of 2-3 weaker than 
in the extreme observed cases can be regarded as a sensible 
approximation. 

As a last application for this geometry we mimic the 
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Figure 12. Same cumulative fraction h^ip as in figure fTTI Here for the metal lines. Positive velocities are infiowing into the galaxy and 
negative velocities are out of the galaxy. Contour lines are at 0.1, 0.01, 0.001 and 0.0001 respectively. Note the different y-axis scalings 
with respect to figure [TTI as well as from row to row. In C ll one sees an infiow > 150 km s"-*^ with an EW > 0.17 A in 0.4% of all 
observations. For Mgll one sees an infiow > 150 km s"-*^ with an EW > 0.2 A in 1.3% of all observations. The line which has by far the 
strongest signal is the Mgll line followed by the Fell line. These values should be achievable by future observations. 
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SIO stacking procedure (see their figures 6 and 10) they 
used to increase the signal to noise ratio. They stacked 89 
galaxies with z = 2.3 ± 0.3 to investigate the kinematics of 
the galaxy-scale outflows. We produce stacked spectra using 
our simulations by summing up the line profiles of several 
thousand different directions for each of the three galaxies 
and stack them together. We determine the absorption line 
profile for a spherical shell between an outer radius and an 
inner radius ri. The outer radius is always kept constant at 
1.0 Rv which corresponds roughly to 74kpc. The inner ra- 
dius n however is varied between 0.3 i?v and 0.02 R^. This 
is done for the following reason: As one can see from fig- 
ure [U there are very clean, unperturbed streams into the 
galaxy down to as far as ^ 0.25i?v In the immediate vicin- 
ity of the galaxy there is the "messy region" without well 
defined streams and with no well defined gas disk either. 
The galactic disk itself has no sharp edge and therefore its 
extent is not determined easily either. In order to account 
for these uncertainties we show several different absorption 
line profiles for different inner radii for every wavelength in 
question. All are averaged in the aforementioned way. In fig- 
ures [13] and [14] the resulting profiles are shown for all the 
lines. These figures can be directly compared to figures 6 or 
10 in SIO. The plots reveal that we have for Lya a FWHM 
?7 of up to ~ 1900 km s~^ with a line depth of up to 0.90 
for an inner radius n going as deep as 0.02 R^ (1.5 kpc). 
The Lya lines for high n having a shallow total MLD show 
very extended wings at the edges of the line profile whereas 
lines for small ri having a deep total MLD the edges of the 
profiles are much sharper since the centre of the line is much 
deeper. All profiles of all lines always peak in the positive, 
indicating inflow. Our strongest metal line is Sill. It has a 
maximum line depth of ^ 0.23 with a rj ^ 250kms~^ for 
ri = 0.02 i?v This confirms our initial guess that metal lines 
are so much weaker than the Lya line, since the inflowing 
material is mainly unprocessed primordial gas with very low 
metallicity. The predicted metal line absorption profiles ap- 
pear tiny compared to the corresponding lines presented in 
figures 6 or 10 of SIO having a line depth of 0.5 and 77 ^ 1000 
km s~^ in metals. Detailed comparisons between our values 
and the observed values by SIO are given in Table [2] 

It is hard to compare our Lya absorption line profiles to 
the observations of SIO as they actually show emission line 
profiles. One could still draw interesting conclusions by com- 
paring the underlying optical depths r of the SIO observa- 
tions with our simulations. However, the Lya emission is not 
only an indicator of cold streams. Probably the most suitable 
lines for the purpose of detecting cold streams in absorption 
are C 11 and Mg 11 since they have the strongest signal. Out 
of those Mgii is closer to being observ ed with the needed 
sensitivity and resolution (cf. for example lMateiek Simco3 
I2OI2I ). As SIO point out, the easiest way for a Lya photon 
to reach the observer is to acquire the velocity of the out- 
flowing material on the far side of the galaxy, and to be 
emitted in the observer's direction. In this case the photon 
is redshifted by several hundred km s~^ relative to the bulk 
of the material through which it must pass to reach us. This 
picture would explain qualitatively why the dominant com- 
ponent of Lya emission always appears redshifted relative 
to the galaxy systemic velocity (see e.g. IPettini et al.ll2002l : 
lAdelberger et al.ll2003l ). implying that their results for Lya 
have to be interpreted with caution. Our results for Lya 
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Table 2. Comparison of maximum line depths (MLD), FWHM 77 
and EW Wq in the central geometry between SlO's observations 
(their figure 10) and our predictions from the simulations (the 
'Vj = 0.02 i?v" line in our figures [T3] and [l4]) . 
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Figure 13. Lyo; absorption line profile for a central source geom- 
etry averaged over different viewing angles and all three galaxies. 
We integrated from 1.0 i?v down to different inner radii r{. Posi- 
tive velocities are inflowing into the galaxy and negative velocities 
are out of the galaxy. Note the different scaling of both axes in 
comparison to figure fl4] The profiles always peak in the positive, 
indicating inflow. Lines for high n having a shallow total MLD 
show very extended wings at the edges of the distribution whereas 
lines for small ri having a deep total MLD the edges of the profiles 
are much sharper since the centre of the line is much deeper. 



absorption in the simulations are not subject to this com- 
plication, therefore one should be cautious when comparing 
our Lya results to the SIO data. Metal lines are less prob- 
lematic. 

There are two main conclusions from our analysis of the 
central-source geometry, which are discouraging concerning 
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Figure 14. Same absorption line profiles for a central source as in figure fTsl this time for the metal lines. Positive velocities are infiowing 
into the galaxy and negative velocities are out of the galaxy. Note the different y-axis scaling from row to row and the different overall axes 
scaling with respect to figure fTsI The metal lines are much weaker than the Lya line, since the infiowing material is mainly unprocessed 
primordial gas with very low metallicity. These metal lines also appear tiny compared to the corresponding lines presented by SIO. The 
profiles always peak in the positive, indicating infiow. The lines with the strongest signals are Si ll and Mg ll. A more detailed quantitative 
comparison to the observations is presented in table [2] 
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Figure 15. Map of the projected neutral hydrogen column den- 
sities (A^Hl) for the background geometry of our fiducial galaxy 
at 2; = 2.3. The box side is 138 kpc (physical). The outer circle 
marks the virial radius and the inner circle is at 1/3 Rv. Con- 
tours are shown for 10^^, 10^^ and lO^^cm"^. One can see the 
three major streams outside 1/3 i?v as well as clumps which are 
actually satellites with dark matter haloes. Both of them have 
Nm > 10^'^ cm-2. The background has A^hi < 10^^ cm-^. in- 
side 1/3 Rv there is the messy region. It has column densities 
with values between those of the stream and the background. 

the potential for detecting cold streams. First stacking of 
absorption line profiles does not strengthen but weakens the 
signature of cold streams this might be counter intuitive, but 
compare figure [6] with figures [T3l and [T4|) . See figures [TTI and 
fT2l for the predicted strengths of a cold stream absorption 
signal when observing without stacking. Second, the Lya 
signal detected by observations is naturally dominated by 
outflows which are ex pected to have much higher cov ering 
factor and metallicity (jFaucher-Giguere KeredlioTTI ). 



5 BACKGROUND SOURCE 

In this section we consider the Lya and metal line absorption 
that occurs as UV light emitted by a background galaxy or 
a quasar is absorbed by gas in the circum-galactic environ- 
ment of the foreground galaxy under investigation. In this 
geometry one cannot distinguish between inflowing and out- 
flowing gas. However, the distribution of impact parameters 
allows us to explore the spatial distribution of the circum- 
galactic gas surrounding the foreground galaxies. SIO con- 
structed a sample of 512 close angular pairs, (1 — 15^^), of 
galaxies in the redshift range z ^ 2 — 3 with a large red- 
shift difference such that one is at the background of the 
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Figure 16. Area weighted average column density Nai as a func- 
tion of impact parameter b for H I and five different metal lines 
in the background source geometry. The Hi column density de- 
creases through the interaction region in the greater disc vicinity, 
from 10^^ cm^ near the centre of the galaxy to 10-^° cm^ at 0.5 
Rv. It remains roughly constant at 10-*^^ cm^ in the stream regime 
out to 1.5 Rv, and is not much lower out to 2.5 Rv 



other with no physical association. The pair separations cor- 
respond to galactocentric impact parameters in the range 
3 — 125 kpc (physical) at z = 2.2, providing a map of 
cool gas as a function of galactocentric distance for a well- 
characterised population of galaxies. The discussion in this 
section will lead to a prediction of absorption line profiles 
that are directly comparable to the ones observed and pub- 
lished in section 6 of SIO. 

For a first impression of the simulation data in the 
background-source geometry we show in figure [15] a map 
of projected neutral hydrogen column density. The box of 
140 kpc is centred on our fiducial galaxy of Mv ^ 4x lO^^M© 
at z = 2.3, which serves as the foreground galaxy. The col- 
umn density NAi{x,y) of element A and ionisation state i 
at the position {x,y) itself is now calculated by 

/bs 
XAi{x) nA{x) dz, (7) 
-bs 

where bs indicates the box size of twice the virial radius, 
nA(^) is the total gas density of element A at position 
X — {x,y,z) and xai is the ionisation fraction of element 
A in state i. This equation is the analog of equation [1] for 
the central geometry. The map in figure [15] corresponds to 
figure [3] in the central geometry. The outer circle marks the 
virial radius and the inner circle is at 1/3 Rv Contours are 
shown for 10^^ 10^^ and 10^^ cm . One can see the three 
major streams which are very pronounced outside 1/3 Rv 
They have neutral hydrogen column densities of A^hi > 10^^ 
cm~^. The background outside 1/3 Rv has neutral hydrogen 
column densities of A^hi < 10^^ cm~^. In the messy region 
inside 1/3 i^v one cannot distinguish between streams, galac- 
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Figure 17. Area covering fraction gc as a function of impact parameter b and cumulative column density for different ionisation states i 
of the elements A in the background source geometry. The panels show for every impact parameter bin b the fraction of the area within 
this impact parameter bin which has a column density of A^Ai or above. The contours mark area covering fractions gc of 1.0, 0.1 and 
0.01. The values are averaged over three different galaxies and over three orthogonal directions. We see that Hi shows the highest column 
densities over the whole radius range. It also reaches a area covering fraction of 1.0 at lO-*^"^ cm~^ over the whole radius range. The 0.1 
and 0.01 area covering fraction lines a decreasing monotonically with radius for all panels. The metal i onisation states with the highest 
area covering fractions are Cll and CiV. T he Hi pan el agrees with the results of Fll (their figure 4), iFaucher-Giguere Kerej (|201lL 
their Table 1) as well as with the results of I S hen et al., (^2012. ) (their figure 10). 



tic disk or background. The neutral hydrogen column den- 
sities within this region vary between 10^^ cm~^ and 10^^ 
cm~^. One should note that this ideal map cannot be di- 
rectly compared to any observed galaxy, because it requires 
a background point source along the same line of sight of 
each pixel in the map. 

The impression from figure [15] is that the column den- 
sity tends to decrease with increasing distance from the 
galaxy centre. Figure [16] shows the area- weighted average 
column density, for Hi and five different metal lines, as a 
function of the impact parameter b. The H I column density 



decreases through the interaction region in the greater disc 
vicinity, from 10^^ cm~^ near the centre of the galaxy to 10^° 
cm~^ at 0.5 Rv It remains roughly constant at 10^^ cm~^ in 
the stream regime out to 1.5 i?v, and is not much lower out 
to 2.5 Rv The metal ionisation states also decrease inside 
0.5 Rv and are roughly constant at values between 10^^ cm~^ 
(Oi) and 10^^ cm~^ (Siiv) in the outer halo and beyond. 

Figure [15] is limited to the angle- averaged column den- 
sity. Figure [T71 presents the distribution of column densities, 
for the purpose of evaluating the likelihood of an observa- 
tion through a given line of sight. We plot area covering 
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Figure 18. Example Lya absorption line profiles seen from a single direction in the background geometry integrated from 0.3 to 1.0 R^. 
For a background source one cannot tell from the redshift or blueshift of the line relative to the central restframe whether it represents 
infiow or outfiow with respect to the galaxy centre. From inspection of the three-dimensional information in the simulations, we know 
that almost in every case a significant coherent shift of the line is associated with infiow. Values for neutral hydrogen column density 
and EW are quoted in each panel. Uppermost panel: various modes of bimodal streaming, intermediate upper panel: different examples 
of strong unidirectional streaming, intermediate lower panel: low absorption profiles and lowermost panel: extraordinary cases. Some of 
the lines are saturated like B, C, D, E, F, K and L, whereas others are not. The letters correspond to the indicated position of the map 
from figure [201 
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Figure 19. Same figure as figure [TS] panel K (the panel with the deepest signal) this time for selected metal lines providing also the 
deepest signal, namely: C ll (1334 A), Sill (1260 A), Mgll (2796 A) and Fell (2383 A). Panel K is fully saturated in Lya but still shows 
peak abs orption line d e pths of 0.25 and EW of ~ 0.2 A for these metal lines. The C ll panel of this plot shows slight differences with the 
results of iKimm et al.l (|201ll , their figure 2, lower panel) due to different line computing algorithms. 



fractions Qc as a function of impact parameter b and cu- 
mulative column density A^Ai- The values are averaged over 
three different galaxies and three orthogonal projections per 
galaxy. The panels show for every impact parameter b the 
fraction of the area within an annulus of radius b and width 
Sb which has a column density of Nai or above. One can see 
that there are obviously much higher columns of H i present 
than of the other metal ionisation state. H i has column den- 
sities which are several orders of magnitude higher, depend- 
ing on the actual radius and the area covering fraction one 
wants to look at. The columns of the metal ionisation states 
lie closely to each other and decrease rather smoothly with 
increasing radius. The = 0.1 lines of the two carbon lines 
are fairly similar, both fall off almost linearly from ^ 10^^ 
cm~^ at the left edge of the graph to ^ 10^° cm~^ at the 
right edge of the graph. This is also true for the Qc — 0.1 
lines of the two silicon lines. They fall off from around 10^^ 
to below 10^. The Qc = 0.01 line, which lies much higher, on 
the other hand has more differences between the low-ionised 
metals (like Cii, Sin or Fell) on the higher ionised metals 
(like C IV or Si iv): In case of the low- ionised metals the line 
has considerably more features like the extended bumps at 
impact parameters of 1.0 or 2.0 R^. The distribution 
for C II indicates the highest, the one for Mgii indicates the 
lowest overall column densities. 

The H I panel of figure \T7\ is comparable to figure 4 of 
Fll, who implemented a full radiative transfer calculation. A 
comparison between the two allows us to evaluate the accu- 
racy of our simplified treatment of self-shielding. For impact 
parameter other than b = 2.0 R^ the latter show in their fig- 
ure 4 values which are normalised to the value it would have 
at b = 2.0 i?v In the following we will only quote the un- 
normalised values. In our results the same covering factors 
correspond to column densities that are usually smaller by 
a factor of ^ 3. For example at 5 = 2.0 R^ Fll have a cover- 
ing factor of 0.01 at 3 x 10^° cm"^ and of 0.1 at 10^'^ cm"^ 
In our results these covering factors correspond to column 
densities of 10^° cm~^ and 3 x 10^^ cm~^, respectively. The 
H I panel of figure pT| is also r oughl y comparable to Table 1 
of iFaucher-Giguere Kered (|201lh . In our results we have 
covering factors which are on average smaller by a factor of 
1.5. For example at b = 0.5 Rv they have an average cover- 




Wo [A] 

Figure 20. Map of the Lyo; absorption line EWs for the back- 
ground geometry of one of our simulated galaxies. The box side 
is 276 kpc (physical). The outer circle marks the virial radius and 
the inner circle is at 1/3 i^v The contour lines indicate an EW of 
1 and 0.1 A respectively. The circles with the letters attached to 
it indicate the position of the example line profiles from figure [TS] 
The three major streams as well as clumps (satellites with dark 
matter haloes) have Wo > 1 A. The background has Wo < 1 A. 



ing factor of - 0.3 at 10^^-^ ^ A^hi ^ 2 x 10^° and of - 0.08 
at > 2 X 10^° cm~^. Our corres ponding c o vering factors are 
0.2 and 0.1. In their figure 10. Shen et all (|2Q12l l quote a Hi 
covering factor of column densities higher than 10^^'^ cm~^ 
for r < i?v (r < 2R^) of 27% (10%). Our corresponding cov- 
ering factors are 0.2 and 0.5. We conclude that our results 
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based on the simplified treatment of self-shielding (section 
O provide a rather good approximation to the results from 
the full radiative transport analysis. 

We next compute example absorption line profiles for 
a background source, in analogy to the line profiles for a 
central source, using equations (|2|) to (|5}. From now on, 
we apply a Gaussian point spread function with a beamsize 
(= FWHM) of 4kpc, following the algorithm described in 
section |4l In addition, the velocity resolution is now down- 
graded to 50 kms~^, to match the velocity resolution used in 
the corresponding figures of SIO. In figure [TSlandllQI we show 
a selection of absorption line profiles, for Lya and metal 
lines, respectively, demonstrating a large variety. These fig- 
ures are the analogs of figures [6] and [3 for a central source. 
Some line profiles show a single dominant stream, others in- 
dicate two streams, some display weak absorption and others 
show strong absorption. For a background source one can- 
not tell from the redshift or blueshift of the line relative to 
the central restframe whether it represents inflow or outflow 
with respect to the galaxy centre. From inspection of the 
three-dimensional information in the simulations, we know 
that almost in every case a significant coherent shift of the 
line is associated with inflow. In the case of Lya, 7 out of 12 
lines in the sample are saturated (B, C, D, E, F, K, L). Panel 
K, which shows the strongest absorption in Lya, has for se- 
lected metal lines a maximum absorption depths of only 0.25 
and an equivalent width EW of only ^ 0.2 A. Some of the 
weaker Lya lines (e.g, panels I or even H) might be com- 
pletely erased by noise and the ISM component ai v — 0. 
The leftmost panel of figure [T9l i s ana logous to the lower 
panel of figure 3 in iKimm et al.l (|201lh . The latter obtain 
for C II a maximum line depth of ^ 0.9 with a FWHM 77 
of 650 km s~^. The difference from our result most likely 
stems from the fact that they used a simulation with a res- 
olution low er by a factor of 20, the Horizon MareNostrum 
simulation (|Ocvirk. Pichon Tevss ier 2008), and addressed 
haloes more massive by a factor of a few. Furthermore, they 
used a different prescription for the Gaussian velocity distri- 
bution, where they obtained their velocity dispersion from 
the neighbouring 26 cells instead of the more appropriate 
algorithm described in equations ([2|) to Q. 

As for the central source, we compute the equivalent 
widths from the absorption line profiles using equation (|6|). 
The EWs for a background source are mapped for our fidu- 
cial galaxy in figure 1201 The outer circle marks the virial 
radius and the inner circle is at 1/3 i?v • This plot is the ana- 
log of figure [8] for a central source. The EWs shown in it 
span two and a half orders of magnitude. The three major 
streams are fairly prominent, in addition to the central re- 
gion in the greater vicinity of the disc. Also very prominent 
are the parts of the IGM which are stream-less. 

Figure [20l illustrates the dependence of the EWs on the 
projected distance from the centre, or impact parameter b. 
In the galaxy vicinity at the halo centre, b < 10 kpc, say, we 
see the highest absorption, which gradually declines toward 
larger radii, as the volume becomes dominated by the dilute 
medium between the streams. Figure [21] shows the average 
EW as a function of impact parameter 6, averaged over the 
three simulated galaxies at z = 2.3, where the virial masses 
are Mv :^ 4 x 10^^ Mq. For Lya, we show for comparison the 
observational results of SIO (their Table 4), assuming a virial 
radius of 74 kpc. In our simulation results for all the lines. 



the inner galaxy is responsible for a steep peak in EWs in the 
range b ^ 0.5i?v Outside of this radius the EW is lower and 
it shows an overall gradually declining with b. A comparison 
with the observations shown in figure 21 of SIO indicates that 
the observational EW for the Lya line lies a factor of two 
or less below our predicted value at low b, and is consistent 
with our prediction at larger radi i. The simulation r esults in 
figure 14 of Fll and in figure 9 of IShen et all (I2OI2I ) for Lya 
are in excellent agreement with our current results, and so 
are the Fll predictions for Sin (1260 A). The observational 
results of SIO for the EW of the metal lines are two orders of 
magnitude above our predictions. These are likely to reflect 
massive, cold, dense and metal-rich outflows that are not 
reproduced in an appropriate amplitude in our simulations. 

The probability to observe a given EW along a given 
line of sight of impact parameter b is shown in figures [22l 
and 1231 Shown in colour is the area covering fraction Qc as 
a function of impact parameter b and cumulative EW for 
different absorption lines. The panels show for every impact 
parameter bin b the fraction of the area within this impact 
parameter bin which has an EW of the indicated value Wo 
or above. The contours mark fractions of 1.0, 0.1, 0.01 and 
0.001. The values are averaged over three different galaxies 
and three orthogonal projections, in both directions each. 
This figure is the analog of figure [TOl for a central source. We 
see that in Lya we have EWs > 0.01 A everywhere. Peak 
EWs go as high as 10 A. The 0.1 area covering fraction line 
is falling monotonically with b, whereas the 0.01 and 0.001 
lines show noisy distortions at b > 1.5 i?v For the metal 
lines the peak EWs lie near 0.1 A and their minimum values 
are well below 10~^ A. In all cases the EWs are declining 
mainly monotonic with radius. For Lya the decline is rather 
gradual, only by an order of magnitude over the studied 
radius range. For the metal lines the decline is by more than 
five orders of magnitude. We see some noisy distortions in 
the high 6, low Qc regime. Over all the strongest lines are 
C II and Mgii. 

Finally, we try to produce absorption line profiles that 
could be compared to those presented in figures 6 and 10 
of SIO. They stacked spectra of several galaxies in order to 
increase the signal to noise ratio. Their sample contains 512 
close (1^^ — 15^^) angular pairs of 2; 2 — 3 galaxies with red- 
shift differences that indicate no physical association. The 
sample explores the cool gas as a function of galactocentric 
impact parameter in the range 3 — 125 kpc (physical) at 
z — 2.2. In figures [231 and [25] we show the resulting absorp- 
tion line profiles for Lya and metal lines in our simulations, 
stacked over six principal different viewing directions (three 
orthogonal axis in both directions) . We show the line profiles 
for the same three different impact parameter bins as in SIO: 
10 - 41 kpc, 41 - 82 kpc and 82 - 125 kpc. (Note the differ- 
ent scales in figures [2411251 and from row to row within figure 
[25]) . The lines for 10 — 41 kpc in those figures correspond to 
figure 18 in SIO, the lines for 41 — 82 kpc correspond to 
figure 19 in SIO, and the lines for 82 — 125 kpc correspond 
to figure 20 in SIO. For Lya we have a FWHM (= ry) of 
up to ~ 550 km s~^ and a relative intensity of 40-50% at 
the line centre of the innermost radius bin. It is not easy 
to compare our Lya absorption line profiles to the observa- 
tions of SIO as they actually show emission line profiles. One 
could still draw interesting conclusions by comparing the un- 
derlying optical depths r of the SIO observations with our 
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Figure 21. Average EW Wq as a function of impact parameter 6, averaged over the three simulated galaxies. For Lya, we show for 
comparison the observational results of SIO (their Table 4). The inner galaxy is responsible for a steep peak in EWs in the range 
b ^ 0.5i?v Outside of this radius the EW is lower and it shows an overall gradually declining with b. As expected there is a strong 
negative correlation between impact parameter b and EW. This behaviour is the same for all ten lines. The Lya an d the Sill (1260 A) 
lines agree with the theoretical results of Fll (their figure 14). The Lya line also agrees with the theoretical results of lShen et ahl (l2012l ) 
(their figure 9). The observational results of SIO for hya (their Table 4 and their figure 21) are also in fairly good agreement. 




0.5 1 1.5 2 2.5 

b[Rv] 

Figure 22. Area covering fraction gc as a function of impact pa- 
rameter b and cumulative EW Wq for the Lyo; absorption line 
in the background source geometry. Shown is for every impact 
parameter bin b the fraction of the area within this impact pa- 
rameter bin which has a EW of Wq or above. The contours mark 
fractions of 1.0, 0.1, 0.01 and 0.001. The values are averaged over 
three different galaxies and the six principal axis. We see that we 
have an EW > 0.02 A everywhere. The 0.1 area covering fraction 
line is falling monotonically with 6, whereas the 0.01 and 0.001 
lines show noisy distortions at b > 1.5 i?v 



simulations. However, due to the special transfer behaviour 
of Lja mentioned at the end of section [4l the Lya emission 
is not only an indicator of cold streams. Therefore the re- 
sults have to be taken with caution. Our strongest metal 
absorption lines. Sin and Cii, have both a maximum line 
depth of - 0.015 with ry - 350 km s"^ for the 10 - 41 kpc 
bin. This appears tiny compared to the corresponding lines 
presented in figures 18 or 19 of SIO having a line depth of 
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Table 3. Comparison of maximum line depths (MLD), FWHM 
T] and EW Wo in the background geometry between SlO's obser- 
vations (their Table 4, second line) and our predictions from the 
simulations (the "10 — 41 kpc" line in our figures [24l and [25)) . 

0.2 and ry - 1000 kms"\ Detailed comparisons between our 
values and the observed values by SIO are given in Table O 
Figure [24] as well as the upper middle panel of figure 
1251 can be compared to figure 13 of Fll. For Lya, in the 
10 — 41 kpc bin, the maximum line depth is 0.61 and 0.55 
respectively, and the FWHM 77 is 550 and 400 km s~^ re- 
spectively. For Sin (1260 A), in the same 10 — 41 kpc bin, the 
maximum depth is 0.013 and 0.05 respectively, and 77 ^ 350 
km s~^ in both. These results are similar. The upper left 
panel of fi gure [25l can be co mpared to the right panel of 
figure 3 in lKimm et al.l (|201lh , where for C 11 the maximum 
line depth is - 0.65 with a FWHM ry of 550 km s"\ We 
identified th r ee po ssible reasons for these differences: First 
iKimm et al.l (|201lh used a simulatiorQ having a much lower 

^ the Horizon MareNostrum simulation 
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b [RJ b [RJ 

Figure 23. Same area covering fraction gc as in figure [22] this time for the metal fine absorption fines. Peak EWs go as high as 0.1 A 
and the minimal values are well below 10~^ A. In all cases the EWs are decreasing mainly monotonic. The decrease goes over five orders 
of magnitude and more. We see some noisy distortions in the high 6, low Qc regime. Over all the strongest lines are C ll and Mg ll. 
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Figure 24. Absorption line profiles for the background geome- 
try for hya. The data are averaged over three different impact 
parameter bins, three different galaxies and six different view- 
ing directions in the background geometry. Note the different x- 
and y-scales in this figure in comparison to Figure [25] This fig- 
ure agrees with the results of Fll (their figure 13). A qualitative 
comparison to the observation is presented in table O 



resolution than we do. Second, the haloes they analysed are 
in a very different (much higher) mass range (Mv > lO"*^^ 
M0) than the haloes we use. Third and most importantly 
they use a different prescription for the Gaussian velocity 
distribution as described earlier. We conclude that stacking 
of the absorption lines from a background source tends to 
wash out the cold filament absorption signal (compare fig- 
ure [18] with figures [24] and [25]) . Unlike the case of a central 
source, in the case of a background source the absorption by 
inflowing gas is not well separated from the absorption by 
outflowing gas. 



6 CONCLUSIONS 

Theory, including hydrodynamical cosmological simula- 
tions, tells us that massive galaxies at high redshifts 
were fed by cold gas streams, inflowing into dark- matter 
halo es at high rates along the fllaments of the cosmic 
web (iKeres et al.ll2005l:lDekel k Birnboimll20Q6l : IPekel et al.1 
l2Q09l : iDanovich et al.1 1201^)" but observing this process is 
not straightforward. Here we addressed the expected ab- 
sorption signature from these cold streams. We "observed" 
simulated galaxies for Lya absorption as well as metal lines 
from low and medium ionisation ions . We focused on th e 
absorption line proflles as observed bv lSteidel et al 
using their set of ten different lines and mimicking their 
way of stacking data from several galaxies, for sources that 



(|Ocvirk, Pichon Tevssieil I2OO8I ) with a maximum resolu- 
tion of 1 kpc 



are either in the background or at the centre of the ab- 
sorbing halo itself. The simulations used are zoom-in cos- 
mol o^ical simulations with a maximum resolution of 35-70 
pc (jCeverino, Dekel BournaudI |201QI ). Self-shielding was 
accounted for by a simple density criterion. 

We showed a variety of example absorption line pro- 
flles, as well as averaged absorption line proflles, both for a 
central source and a background source, which are directly 
comparable to the observations of SIO. We computed the 
resulting EWs and their sky covering f ractions for the two 
geometries. We found in agreement with lKimm et all (| 20 111 ) 
that when observations are averaged the absorption signa- 
tures of the low-ionisation metals from the cold streams are 
overwhelmed by the noise of other absorption processes in 
the host galaxy. This has the following reasons: First, since 
the infalling gas has compared to the ISM of the host galaxy 
so low densities and so low metallicities that subsequently 
compared to other absorption processes in the host galaxy 
the low-ionisation transitions of the cold streams also have 
tiny optical depths. Second, the difference in redshift be- 
tween the absorption signal of the cold streams and that of 
the ISM is too small. The general absorption signature of the 
ISM still overwhelms the one from the cold streams even in 
the very redshifted part of the line. Last but not least, the 
flows have a radial, narrow stream geometry, which makes 
them only very occasionally aligned with the line of sight for 
a high optical depth. To show this we computed the proba- 
bilities of detecting an inflow of a given speed with a given 
signal. 

Similar to our work, Fll studied the absorption char- 
acteristics of the gas in galaxies and streams, in order to 
compare with the statistics of observed absorption-line sys- 
tems. Like us they postprocessed high-resolution simulated 
galaxies for determination of the ionisation states of the 
gas. To compute those however we used the combined ef- 
fects of electron-impact collisional ionisation and photoioni- 
sation together with a distinction between UV-shielded gas 
(where collision-ionisation equilibrium was assumed) and 
UV-unshielded gas (where photo-ionisation equilibrium was 
assumed). Fll on the other hand used a full radiative trans- 
fer approach. Since both analyses use the same data, the 
high resolution CDB simulations, direct comparisons to Fll 
gives us a good handle to estimate the strengths and weak- 
nesses of our simplifying procedure. We compared the dis- 
tributions by volume of the neutral hydrogen fractions xm 
and total hydrogen gas densities nu in the circumgalactic 
environment (our flgure[2]). We also compared the Hi area 
covering fractions (our flgure [TT)) as well as the average EW 
as a function of impact parameter b for Lya and the Sill 
(1260 A) lines (our flgure [2T]) both for the background ge- 
ometry. We flnally compared the absorption line proflles for 
the background geometry in Lya and Sill (1260 A) (our flg- 
ures 1241 and I25p . The differences from our results are small, 
justifying the assu mptions made in our simpli f ying m odel. 

Both Fll and iFaucher-Giguere KeresI (l201lh anal- 
ysed the sky covering fractions of neutral hydrogen col- 
umn densiti e s, and came up with similar results to ours. 
iKimm et al.l (|201lh predicted C 11 absorption line proflles 
similar to ou r s. Cer tain differences emerge from the fact that 
iKimm et all (|201lh deploy a Gaussian proflle instead of a 
Voigt proflle. Also, while we used thermal Doppler broaden- 
ing that depends on the temperature of the cell, they used 
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Figure 25. Same absorption line profiles for the background geometry as in figure [24] this time for the metal lines. Note the changing 
y-scales from row to row as well as the different x- and y-scales with respect to figures 1241 T he Sill (1260A) pa nel of this plot agrees 
with the results of Fll (their figure 13) the C ll panel of this plot disagrees with the results of iKimm et al.l (|201ll . their figure 3) due to 
different line computing algorithms. The lines with the deepest signals are Cii and Sill (1260 A). A more detailed qualitative comparison 
to the observation is presented in table O 
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turbulent Doppler broadening that depends on the veloc- 
ity differences with respect to the 26 neighbouring cells. 
Fll, like us, concluded that cold streams are unlikely to 
produce the large EWs of low-ion metal absorption around 
massive galaxies. If the average equivalent widths reported 
by SIO are confirmed (e.g. with higher resolution observa- 
tions), more precise feedback models must be included in 
the simulations, which are capable of causing winds with 
the right densities and temperature range, since only such 
winds are able to produce these high equivalent widths of 
low-ion metal absorption lines. 

We mention three potential limitations of our analysis. 
A limitation of our simulations may arise from the artificial 
pressure floor imposed in order to properly resolve the Jeans 
mass. This may have an effect on the temperature and den- 
sity of the very dense and cold parts of the streams, with 
potential implications on the computed absorption. Still, the 
AMR code is the best available tool for recovering the stream 
properties. With 35-70 pc resolution, and with proper cool- 
ing below lO^K, these simulations provide the most reliable 
description of the cold streams so far. 

SIO and others observed strong signatures of metals at 
distances as large as 200 kpc from the galaxy centre, that are 
interpreted as massive outflows. Such massive and cold out- 
flows are not reproduced in our current simulations, which 
incorporate supernova feedback, but not yet strong radia- 
tive feedback and other processes that are capable of driv- 
ing sufficiently massive outflows. Our studies of absorption 
from the simulated inflows are valid only if one could as- 
sume that the outflows do not drastically affect the inflowing 
streams. The stro ng outflows produced in the simulations of 
IShen et"al ] (120121 ) indeed avoid the incoming dense narrow 
streams, and find their way out in a wide solid angle through 
the dilute medium between the streams. The same is true 
for the weaker and hot outflows that are produced in our 
simulations, with a mass loading factor of about a third. 
Therefore our analysis is reliable. 

The ionisation states should ideally be coupled to 
the hydrodynamical calculations using radiative transfer, 
but in our current analysis they are computed in post- 
processing. The analysis distinguished between UV-shielded 
gas, where collision- ionisation equilibrium was assumed, and 
UV-unshielded gas, where photo-ionisation equilibrium was 
assumed. An additional source of local ionising radiation is 
the Lyman continuum from newly born stars. This could be 
a non-negligible source of ionisation, which is neglected in 
our calculations so far. 

We conclude that the signatures of cold inflows are sub- 
tle, and when stacked are overwhelmed by the outflow sig- 
natures. Our predicted Lya line absorption profiles agree 
with the observations, while the stacked metal line absorp- 
tion from the inflows is much weaker than observed in the 
outflows. The single-galaxy line profiles predicted here will 
serve to compare to single-galaxy observations. 
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